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ESTIMATED  PROPERTIES  OF  LAMINAR  WAKES  WITH 
COUPLED  RATE  CHEMISTRY  AND  DIFFUSION 


M.  H.  Bloom,  R.  W.  Byrne  and  E.  D.  Kennedy 

ABSTRACT 

Theoretical  analyses  bearing  on  wake  observables  required  a  careful 
representation  of  the  rate  chemistry  of  air  and  its  coupling  with  diffusive 
processes.  A  selection  of  preliminary  results  of  recent  calculations  of 
laminar,  constant-pressure  axisymmetric  wakes,  including  chemical- 
diffusive  coupling  made  for  several  full-scale  flight  and  firing  range 
conditions  is  presented.  The  analysis  employs  an  integral  method  wherein 
the  governing  differential  equations  are  satisfied  exactly  along  the  axis. 

The  present  chemical  and  thermodynamical  models  which  consider  N2, 

N,  02,  O,  NO,  NO+  and  e”  species  are  more  detailed  than  those  usually 
used  in  connection  with  wake  calculations.  However,  the  following 
assumptions  are  retained:  a.  vibrational  equilibrium,  b.  Fick's  law  of 
diffusion  for  each  species,  c.  a  single  representative  Lewis  number  is 
taken  constant  and  equal  to  1.4,  d.  a  constant  value  of  the  Prandtl  number 


was  taken  as  equal  to  0.71, 


Attention  is  devoted  to  the  following  items: 

a.  Properties  along  the  axis  for  representative  reentry  trajectories. 

b.  Wake  lengths  associated  with  velocity,  temperature,  and  density. 

c.  Scaling  considerations  pertinent  to  a.  and  b. 

d.  An  indication  of  the  sensitivity  of  the  results  to  changes  in  the 
chemical  rate  constants. 

e.  Sensitivity  of  the  results  to  different  assumed  values  of  the  Lewis 
cind  Prandtl  numbers. 

f.  The  relative  importance  of  diffusive  and  chemical  terms  evaluated 
along  the  axis  in  each  equation  of  species  concentration. 

g.  The  influence  of  a  higher  "effective"  viscosity  on  the  wake 


characteristics. 
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I  INTRODUCTION 

Theoretical  analyses  bearing  on  wake  observables  require  a  careful 
representation  of  the  rate  chemistry  of  air  and  its  coupling  with  diffusive 
processes  (  see  Refs,  1  and  2)  , 

A  preliminary  selection  of  results  of  recent  calculations  of  laminar, 
constant  pressure,  axisymmetric  wakes  including  chemical-diffusive  coup¬ 
ling,  is  presented  here.  The  calculation  is  derived  from  an  integral  method 
analysis  wherein  the  governing  differential  equations  are  satisfied  exactly 
along  the  axis.  The  present  chemical  and  thermodynamic  models,  which 
consider  N2,  N,  02,  O,  NO,  NO*  and  e  species,  are  more  detailed  than 
those  previously  described  in  connection  with  the  wake  calculations  of 
Refs,  1  and  2  (the  revisions  are  described  in  Ref,  3).  However,  the 
following  assumptions  are  retained:  (a)  vibrational  equilibrium,  (  b)  Fick's 
law  of  diffusion  for  each  species,  (c)  a  single  representative  Lewis  number 
taken  constant  and  equal  to  1.45,  and  (  d)  constant  Prandtl  number  equal  to 
0,71,  The  calculations  were  made  with  an  IBM  7090  digital  computer. 
Conditions  at  the  beginning  of  the  far  wake  were  obtained  by  means  of  an 
integration  of  the  species  rate  equations,  and  the  one-dimensional  equations 
of  energy  and  momentum  along  the  body  surface  streamline  from  the  stagna¬ 
tion  point.  The  general  method  of  calculation  is  described  more  fully  in 


Ref,  4. 
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Results  are  presented  of  calculations  about  the  Trailblazer  II  vehicle 
for  the  Lincoln  Laboratory  under  the  Reentry  Physics  Program. 

Attention  is  devoted  to  the  following  items: 

a.  Properties  along  the  wake  axis. 

b.  Wake  lengths  associated  with  temperature  and  electron 
density,  i.e.,  distances  required  for  the  indicated  properties  to  decay  to  a 
specified  level  or  degree  of  nonuniformity. 

c.  Scaling  considerations  pertinent  to  (a)  and  (  b)  . 

d.  An  indication  of  the  sensitivity  of  the  results  to  a  number 
of  changes  in  the  rate  constants  (at  one  flight  condition,  assuming 
unchanged  initial  conditions)  . 

e.  Sensitivity  of  the  results  to  different  assumed  values  of  the 
Lewis  and  Prandtl  numbers. 

f.  The  relative  importance  of  diffusive  and  chemical  terms 
evaluated  along  the  axis  in  each  equation  of  species  concentration. 

g.  The  influence  of  a  higher  "effective  viscosity"  on  the  wake 
characteristics , 


These  items  will  be  taken  up  in  sequence. 
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II.  DISTRIBUTIONS  ALONG  THE  AXIS  AND  "LENGTHS" 

Wake  Lengths  of  properties  dominated  by  diffusive  processes  are  of 
2 

order  6  pu/fi,  where  6  is  a  typical  wake  radius  roughly  proportional  to 
i  / 

(  Cj-jA)  andpu/p  is  evaluated  at  a  significant  reference  condition,  say 
a  mean  value  of  properties  throughout  the  wake.  For  preliminary  simplified 
discussion,  the  reference  Reynolds  number  per  foot  pu/p  will  here  be 
evaluated  in  the  free  stream.  Thus  a  scaling  parameter  for  the  axial 
coordinate  or  for  wake  lengths  of  diffusive  processes  is  62  pu/p  or  C^A 
(p  u/p)  .  For  a  given  body,  this  scale  factor  depends  only  upon  pu/p. 
Moreover,  over  much  of  a  typical  reentry  trajectory,  changes  in  p  are  more 
influential  than  changes  in  u/p — hence  diffusive  lengths  may  be  expected  to 
scale  roughly  like  p  and  increase  monotonically  during  descent.  The  axial 
variation  of  temperature,  velocity  and  electron  density  for  Trailblazer  II  are 
shown  in  Figs.  1,  2  and  3  respectively.  Figs.  4,5,  and  6  show  the  same 
parameters  plotted  in  terms  of  a  scaled  length  (the  reference  is  the  200,000 
ft.  case)  in  accord  with  what  has  been  said  above.  Complete  correlation  on 
this  basis  cannot  be  expected,  particularly  for  properties  of  the  flow  which  are 
significantly  influenced  by  chemical  rate  processes.  For  example,  the  pro¬ 
nounced  temperature  peaks  indicated  in  Fig.  1  are  primarily  the  result  of 
coupled  diffusion  and  recombination  of  nitrogen  atoms,  the  temperatures 
being  particularly  sensitive  to  small  differences  between  the  stagnation 
enthalpy  and  the  sum  of  the  kinetic  energy  per  unit  mass  and  the  chemical 


energy  contribution  to  enthalpy.  Furthermore,  improved  correlation  is  obtained 
when  account  is  takenof  the  differences  in  initial  condit  ions  and  of  the  influence  of 
the  elevated  interior  wake  temperatures  on  the  diffusion  time  6  p/p.  For 
example,  it  has  been  shown  analytically  that  in  the  present  analysis  uQ/ue 

-Xf 

is  primarily  a  function  of  the  coordinate  (pelle^c/p)  (pQ/p^  )dx,  where 

0^=  Cj}A/4tt.  An  approximate  correlation  of  rate  chemistry  effects  on  the 
axial  scale  can  be  expected  by  a  normalization  of  the  axial  coordinate  with 
respect  to  pu/p  due  to  the  fact  that  two-body  collision  processes  take  place 
roughly  in  proportion  to  p,  whereas  three-body  collision  processes,  such  as 
recombinations,  take  place  in  proportion  to  p2.  The  point  here  is  that  var  ¬ 
iations  of  u/p  during  a  trajectory  are  not  as  drastic  as  that  of  p  and  that  at 
least  a  part  of  the  p 2  scaling  required  for  recombination  will  be  provided  by 
the  p  scaling.  At  low  altitudes,  where  equilibrium  chemistry  prevails,  the 
scaling  applies  to  the  velocity  and  enthalpy.  The  distribution  of  species  then 
follows  from  the  equilibrium  conditions. 

Examination  of  the  ],.!ots  of  electron  densities  clearly  shows  that  one 
must  specify  a  species  level  before  drawing  a  conclusion  on  the  variation  of 

9 

an  associated  wake  length.  For  example,  at  an  electron  density  level  of  1  0  , 

the  wake  lengths  first  increase  during  descent,  then  peak  in  the  neighborhood 
of  140,000  feet  in  the  present  instance,  and  finally  decrease  as  chemical 

7 

equilibrium  is  approached.  On  the  other  hand,  at  very  low  levels  like  10 
or  106  elect rons / cc ,  the  wake  length  increases  monotonically  during  descent. 
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’i»e  establishment  of  peak  values  of  the  concentrations  is  responsible  for 
‘  •us  type  of  phenomenon,  and,  as  mentioned  previously,  the  scaling  con* 
siderations  do  not  apply  to  the  peak  values. 

The  temperature  and  electron  density  plots  also  show  equilibrium 
values  (  dashed  lines)  .  The  equilibrium  temperatures  were  obtained  from 
a  Mollier  diagram  for  equilibrium  air  corresponding  to  the  pressure  and 
enthalpy  determined  in  the  nonequilibrium  wake  calculations.  The  equilib¬ 
rium  electron  densities  were  calculated  from  the  equilibrium  of  the  reaction 
N  +  O  NO  +  e  at  the  density,  temperature,  and  atom  mass  fractions 
of  the  nonequilibrium  claculation.  As  would  be  expected  on  the  basis  of  the 
nonequilibrium  calculations  employed  to  determine  the  initial  conditions  of 
the  wake,  the  equilibrium  temperatures  are  much  higher  than  the  non¬ 
equilibrium  values  at  first.  Further  downstream,  the  nonequilibrium 
temperatures,  substantially  affected  by  chemical  processes  at  low  altitudes, 
diffusion  at  high  altitudes,  and  coupled  effects  in  between,  drive  toward 
equilibrium  and  indeed  may  overshoot  the  equilibrium  values.  This  produces 
the  peak  temperatures  shown  in  the  nonequilibrium  cases.  The  decrease  in 
the  initial  slope  of  the  temperature  curve  during  descent  shows  the  decreasing 
relative  influence  of  diffusion  with  decreasing  altitude.  At  the  lowest  altitude 
(  100,  000  ft,  )  ,  this  effect  and  the  very  much  reduced  flight  speed  combine  to 
yield  the  monotonic  trend  toward  equilibrium.  Downstream  of  the  region  of 
peak  temperatures,  the  nonequilibrium  temperatures  approach  the  equilibrium 
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alues,  although  differences  on  the  order  of  20%  or  more  may  still  be  in 
evidence  before  the  very  far  wake  is  reached.  In  the  very  far  wake,  the  non¬ 
equilibrium  temperatures  of  course  approach  the  equilibrium  values. 

Consider  now  the  electron  densities.  In  this  discussion,  the  term 
"equilibrium11  in  regard  to  electrons  refers  to  the  N  +  O  <=*  NO  +  e 
reaction  in  the  nonequilibrium  bath  provided  by  the  other  species.  It  is 
seen  that  at  the  intermediate  altitudes,  the  nonequilibrium  electron  densities, 
which  are  initially  higher  than  the  equilibrium  values,  drive  strongly  down 
toward  equilibrium  and  then  follow  the  equilibrium  values  closely,  peaking 
along  with  them.  At  these  altitudes,  the  distributions  are  clearly  insensitive 
to  the  initial  values.  At  the  lowest  altitude  {  100,  000  ft)  ,  there  is  the  same 
drive  down  to  equilibrium.  However,  in  this  case,  the  equilibrium  values 
exhibit  no  peak,  or  at  most  a  weak  spread-out  peak,  due  to  the  lower  energy 
of  the  flight  condition.  Hence  the  nonequilibrium  electron  density  exhibits 
no  peak. 

At  the  high  altitudes,  the  equilibrium  distributions  do  peak,  but  the 
nonequilibrium  values  which  are  strongly  influenced  by  diffusion  and  drive 
down  toward  the  equilibrium  values,  may  or  may  not  intersect  the  equilibrium 
peaks.  The  reason  for  this  is  that  the  nonequilibrium  temperatures  to  which 
the  "equilibrium"  electron  densities  correspond,  are  simply  not  high  enough 
to  produce  peak  equilibrium  electron  densities  of  sufficient  magnitude  to 
meet  the  initially-high  nonequilibrium  values.  As  a  result  of  this  (  as  well  as 
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6  'ong  diffusion)  peaks  in  nonequilibrium  electron  densities  are  either  small 
or  nonexistent  at  high  altitudes* 

Another  indication  of  the  possibility  of  scaling  the  wake  results  is 
given  in  the  comparison  betwe’en  the  full-scale  Trailblazer  II  wake  results 
calculated  for  the  150,000-foot  altitude  case  with  the  wake  results  obtained 
for  a  pellet  of  the  same  geometry  but  reduced  in  size  by  a  factor  of  24,5  and 
fired  in  a  range  at  approximately  the  same  ambient  conditions.  The  pellet 
firing  range  conditions  for  which  the  wake  was  calculated  correspond  to  a 
velocity  of  22,  100  fps  and  a  pressure  of  1  mm  Hg*  These  correspond  closely 
to  the  full-scale  conditions  for  the  150,  000-foot  altitude  case  for  the  full- 
scale  Trailblazer  II  (  Table  I)  .  A  comparison  of  electron  densities  along 
the  wake  axis  is  shown  in  Fig*  7.  Major  diffusion  effects  are  extracted  by 
use  of  the  coordinate  scale  (  Di  /  D2)  2x  for  the  small  pellet.  Since  nonequili¬ 
brium  effects  over  the  body  have  been  taken  into  account  in  establishing  the 
initial  conditions  in  both  cases,  the  initial  values  shown  in  Fig*  7  differ  by 
an  order  of  magnitude,  the  more  frozen  pellet  values  being  higher  than  those 
of  the  larger  body.  To  remove  some  effect  of  the  difference  in  initial  con¬ 
ditions,  the  pellet  values  are  also  shown  divided  by  10  (  dashed  curve)  r  The 
short  diffusion  times  in  the  pellet  wake  prevent  the  formation  of  a  peak  in 
electron  density,  despite  the  fact  that  significant  temperature  rises  are 
produced  simultaneously  due  to  diffusion  of  the  major  atomic  species.  The 
influence  of  chemical  effects  is  also  exhibited  by  the  fact  that  the  electron 
densities  of  the  large  body  are  lower  than  those  of  the  pellet.  This  is  so 
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eyond  the  peak  region  even  when  the  influence  of  the  initial  value  s  is 
normalized  out  (dashed  curve)  .  This  trend  is  consistent  with  that  shown 
in  Fig.  3  which  shows  equilibrium  electron  densities  (  at  the  nonequilibrium 
stat^  lower  than  those  given  by  the  nonequilibrium  electron-producing  reaction. 


III.  SENSITIVITY  TO  RATE  CONSTANTS 


Within  the  framework  of  detailed  calculations  of  wake  chemistry  it  is 
of  interest  to  inspect  the  influence  of  changes  in  values  of  the  rate  constants, 
since  the  values  used  originally  and  incorporated  in  the  IBM  7090  program 
were  those  suggested  in  Ref.  5  as  being  within  the  current  state  of  the  art. 
They  are  admittedly  crude  estimates  in  some  cases  and  are  subject  to 
further  revision.  Newly  suggested  values  (  Ref.  6)  differ  from  the  former 
values  in  the  temperature  range  of  interest,  in  some  instances  by  factors 
on  the  order  of  5.  In  the  present  calculations  a  number  of  changes  in  the 
rate  constants  were  made  using  the  important  150  Kft  flight  condition  for  the 
Trailblazer  II  trajectory. 

Reaction  rates  for  the  following  three  reactions  were  varied  by  factors 
of  10  and  10  1  from  the  original  values  found  in  Reference  5. 

a)  N2  +  X  +  9«  8  ev  N  +  N  +  X  (  N2  catalyst) 

b)  N  +  O  +  2.  8  ev  NO+  +  e 

c)  02  +  X  +  5.  1  ev  O  +  O  +  X  (  Ar,  N,  N^,  NO  catalysts) 

The  corresponding  reaction  rates  from  Ref.  5  and  used  in  Ref.  3  are 

5  ,  -1.5  6  2 

Kr  =  .  5  x  10  x  (  T/4500)  ,  cm  /mole  sec 

a 

-3  -3/2  3, 

Kr  ^  =3x10  x  T  ,  cm  /sec/part 

Krc=  2/3  x  10^  x  (T/4500)  ^  cm^/mole2sec 

The  rates  Kr,,  Kr,  ,  Kr  were  varied  one  at  a  time  to  give  the  results 
a  b  c 

shown  in  Figure  8,  where  the  axis  electron  density  values  in  the  wake  are 


given  for  each  of  the  cases  shown  in  Table  II. 


It  is  emphasized  that  the  initial  conditions  which  would  also  be 
subject  to  change  due  to  changes  in  reaction  rates  were  kept  the  same  for 
the  present  test.  It  is  seen  that  the  changes  in  electron  densities  in  general 

are  not  drastic.  The  only  significant  change  in  the  results  is  brought  about  by 

_  \ 

slowing  down  the  electron-producing  reaction  by  a  factor  of  10  (  Case  5)  . 

Increases  in  the  rate  constants  are  less  influential  than  this  decrease  in  the 


electron-forming  rate  constant. 
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IV.  LEWIS  AND  PRANDTL  NUMBER  EFFECTS 

Indications  of  the  importance  of  the  values  of  the  Lewis  and  Prandtl 
numbers  on  the  diffusion  of  species  and  on  the  temperatures  have  been  given 
in  References  1  and  2.  The  sensitivity  of  the  temperatures  to  these  factors 
is  due  primarily  to  the  small  enthalpy  differences  cited  previously,  which 
are  affected  significantly  by  species  concentrations  through  their  effect  on 
the  reaction  energy. 

Wake  calculations  were  made  for  three  sets  of  values  of  Lewis  and 
Prandtl  number  for  the  150,000  ft  altitude  case  on  the  Trailblazer  II  trajectory. 
The  values  of  the  Lewis  and  Prandtl  number  used  throughout  the  calculations 
discussed  above  were  1.45  and  .71  respectively.  The  two  additional  pairs 
of  values  for  which  the  wakes  were  calculated  are  (  Le  =  1.0,  Pr  =  0.71) 

and  (  Le  =  1 . 0,  Pr  =  1.0). 

The  results  of  these  calculations,  plotted  in  Figure  9  show  that  the 
temperature  peak  is  lowest  when  the  Prandtl  and  Lewis  numbers  are  unity. 

The  peak  is  raised  when  the  Prandtl  number  is  0.71  and  the  Lewis  number 
is  unity.  The  highest  temperatures  are  obtained  when  the  Prandtl  number 
is  0.71  and  the  Lewis  number  1.4.  The  corresponding  electron  densities 
given  in  Figure  10  exhibit  the  same  trend. 
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V.  RATIO  OF  DIFFUSION  TO  CHEMICAL  PROqUCTION  TERMS 

To  assist  in  understanding  the  indicated  trends  more  deeply  the 
ratio  of  the  diffusion  to  the  chemical-rate  term  in  each  species -conservation 
differential  equation  was  calculated.  These  are  shown  in  Figure  11  for  the 
representative  Trailblazer  II,  150  Kft  case  for  the  initial  portion  of  the  wake. 

In  this  figure  the  dashed  line  indicates  that  the  diffusive  term  is  negative  for 
the  given  species,  hence  the  sign  of  the  chemical  term  can  be  deduced  from  the 
figure. 

An  interesting  result  is  that  the  chemical  effect  on  the  electron-producing 
reaction  is  dominant  at  the  start  of  the  wake  and  probably  diminishes  in  impor¬ 
tance  downstream.  It  is  reasonable  to  expect  that  diffusion  will  dominate  the 
wake  behavior  over  the  large  downstream  extent. 

More  detailed  studies  of  the  type  indicated  are  in  preparation. 
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rL,  EFFECT  OF  ASSUMING  AN  AUGMENTED  MAGNITUDE  OF  VISCOSITY 

In  studies  of  the  effect  of  turbulence  on  mean  propertie s ,  m  any 
arbitarary  assumptions  must  be  made,  particularly  with  regard  to  turbulent 
transports  in  flows  with  thermal  and  chemical  effects.  In  this  spirit  and 
because  a  calculation  would  be  made  conveniently,  it  was  decided  to  simulate 
very  crudely  the  influence  of  turbulence  in  the  coupled  diffusive -chemical- 
rate  calculation  by  simply  altering  the  magnitude  of  the  laminar  molecular 
viscosity  coefficient.  The  diffusive  and  thermal  transpprts  were  in  effect 
augmented  in  the  same  proportion  since  the  Prandtl  and  Lewis  numbers  were 
kept  the  same  as  the  laminar  values.  Thus  an  ''effective  viscosity”  independent 
of  velocity  gradient  was  introduced.  A  set  of  results  for  electron-density 
distributions  according  to  this  rough  model  is  shown  in  Figure  12,  where 
the  laminar  (Ji  is  multiplied  by  factors  of  10  and  100. 

The  effect  of  this  assumption  is  to  decrease  the  diffusion  time  and 
Reynolds  number  while  the  reaction  times  and  particle  residence  times  are 
unchanged.  This  would  not  be  the  case  if  factors  like  density,  velocity  or  body 
size  were  altered.  The  results  show  the  sharply  reduced  electron  levels.  This 
reduction  is  caused  not  only  by  simple  diffusion,  but  also  by  the  coupled 
reactions  as  is  demonstrated  by  the  scaled  representations  shown  dashed. 

The  scales  of  the  dashed  curves  are  indicated  multiples  of  the  physical 
scale,  and  would  correlate  a  purely  diffusive  phenomenon. 
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TABLE  I 

FLIGHT  CONDITIONS 


Trailblazer  II 

Altitude 

(  thousands  of  ft) 

100 

150 

200 

250 


Velocity 
(  fps) 

15,  875 

21, 600 

22, 475 

22, 550 


TABLE  II 


VARIATIONS  OF  REACTION  RATES  USED 


1. 

no  c 

hange 

2. 

Kr 

a 

x  10 

3. 

Kr 

a 

x  10’ 

4. 

Krb 

x  10 

5. 

Krb 

x  lo" 

6. 

Kr 

c 

x  10 

7. 

Krc 

x  10' 
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NOTATIONS 
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reference  area  for  drag 
Drag 


body  drag  coefficient  (  D/ A) 
electron  density,  electrons 
streamwise  velocity 
streamwise  coordinate  (  ft) 
wake  radius 
density 

absolute  viscosity 


/  cc 
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Subscript 

c  free-stream  (undisturbed)  conditions 
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axis  conditions 
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MODEL  DATE 


Distance  Behind  Body,  ft. 

Figure  3  .  Variation  of  Axis  Value  of  Electron  Density  with  Distance  Behind  Body 

(  Trailblazer  II) 
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Figure  6  Variation  of  Axis  Value  of  Electron  Density  with 
Scaled  Distance  Behind  Body 
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MOdEL  DATE 
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Figure  7  „  Comparison  of  Calculated  Wakes  for  Pellet  and  Full-Scale 
Vehicle.  Axis  Values  of  Electron  Density  Versus  Distance 
Behind  the  Body 
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Fig.  12 


